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AH APPLICATION OE THE VON KARMAN- MILLI&AN 
LAMINAE- BOUNDARY- LATER THEORY AND 
COMPARISON WITH EXPERIMENT 
By Albert E. von Doenhoff 



SUMMARY 



The von Karman-Killlkan theory of laminar "boundary 
layers presented in 1 T .A.C.A. Technical Report No* 504 is 
applied to the laminar "boundary layer about an elliptic 
cylinder on which "boundary- lay er and pressure-distribution 
measurements were made at the National Bureau of Standards. 
An outline of the procedure of the von Karman-Millikan 
method is given. 

Good agreement is obtained between the calculated and 
experimental results, indicating that the method may be 
applied generally to the laminar boundary layer about any 
body provided that an experimentally determined pressure 
distribution is available. It appears that for all 
Reynolds Numbers above 24,000 the separation point for the 
elliptic cylinder should occur at a constant distance be- 
hind the point of minimum pressure, provided that the 
boundary layer does not become turbulent. 



INTRODUCTION 



From consideration of the relative order of magnitude 
of the various terms in the fundamental differential equa- 
tions of flow, Frandtl has shown that the effects of vis- 
cosity at the high Reynolds Numbers common to aeronautical 
problems are of appreciable importance only in the thin 
fluid layer next to the surface of the body, thus consid- 
erably simplifying the fundamental equations. The actual 
solution of the boundary- lay er equations has been a matter 
of considerable difficulty. It was shown, however, that 
separation of the flow from the surface is to be expected 
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when the boundary layer encounters rising pressures in the 
direction of flow. 

The simplest case, that of the boundary layer along a 
flat plate with zero pressure gradient, has "been success- 
fully treated by Blasius. For this case no separation was 
to he expected. The treatment of the flow, about bodies 
when pressure gradients exist has thus far not "been com- 
pletely successful. Numerous methods have been devised 
but either the prediction of tne location of the separa- 
tion roint has been of doubtful accuracy or tne methods 
have been applicable only to limited types of pressure 
dist ribut ions. 

Experiments on an elliptic cylinder have been con- 
ducted at the National Bureau of Standards to study the 
separation phenomena and, in particular, to check the the- 
ory of Pohlhausen. (See reference 1.) The results, com- 
pared with Pohlhausen 1 s approximate solution of the von 
Ea'rn*a'n integral equation, have demonstrated that the Fohl- 
hausen method could only be depended upon to give reliable 
results when the pressure in the outside potential flow 
was decreasing; with increasing pressure separation may 
actually occur when Pohlhausen ? s method fails to indicate 
any flow separation. 

Von Karma'n and hillikan have recently devised a theo- 
ry of laminar boundary layers involving separation (refer- 
ence 2) that appears to have certain advantages over pre- 
vious theories. This solution of the laminar boundary- 
la;v er problem can be applied more generally than former 
methods; the equation that is used to determine the sepa- 
ration point is a close approximation near the surface 
where the separation condition is applied, thus permitting 
a more accurate det erminat ion of the separation point. 

The purpose of the present study is to investigate 
the extent to wnich the theory presented in reference 2 
may be expected to ^ive reliable results. This investiga- 
tion was accomplished by applying the theory to the bound- 
ary layer about the elliptic cylinder on which the Bureau 
of Standards tests were made. A comparison is made between 
the calculated and experimental results, especially with 
reference to the location of tne separation point. 
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EXPERIMENTAL DATA 

The tests (reference l) consisted of pres sure-di st r i- 
bution and boundary- lay ex measurements on an elliptic cyl- 
inder. The pressure dist ribttt ion was measured "by a manom- 
eter connected to orifices in the surface of the cylinder, 
A hoc-wire anemometer was used to make the boundary- lay er 
surveys. The elliptic cylinder on which the tests were 
made had a major axis of 11.78 inches and a minor axis of 
3.98 inches. It was placed in the air stream with the ma- 
jor axis parallel to the general flow. The tests were 
made at a Reynolds Number of approximately 24,000, "based 
on the length of the minor axis. „ 

OUTLINE OF THE VOU KAEMAN- MI L L IEAN PROCEDURE 

The solution of the boundary- lay er problem presented 
in reference 2 is divided into two parte, an outer and an 
inner solution. The outer solution is obtained by trans- 
forming the boundary- layer equation into a form analogous 
to that for the conduction of heat, the solution of which 
is well known. This part of the solution is most accurate 
in the region where the boundary- lay er velocities are near- 
ly equal to those in the outside stream, i.e., the outer 
part of the boundary layer. For reasons peculiar to the 
analysis, it is convenient to join the two solutions at an 
inflection point of the boundary- lay er velocity profiles. 
Thus the range of this solution is taken as the region 
from the outer part of the boundary layer to the inflec- 
tion point, or to the wall if no inflection point exists. 

The inner solution is obtained by transforming the 
boundary- lay er equation by certain approximations into an 
ordinary differential equation, which is most accurate 
near the wall. The inner solution is used only when the 
boundary- lay er velocity profiles show an inflection point. 
Its range is then over the region from the Wall to the in- 
flect ion po int . 

The steps involved in computing the characteristics 
of the laminar boundary layer by this method follow. 

Defini tion, of dimen aionless variahl e s . - In th e fol- 
lowing computations the velocity at an infinite distance 
from the cylinder U 0 is considered to be the unit of ve- 
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locity « The minor axis 1 of the ellipse is used as the 
unit of length. These units of velocity and length lead 
to the formation of the following dimensionless quantities: 

R ' ss " U -^- , Reynolds ITu.mfce.r (^., kinematic viscosity). 

, • • . Cp y« 

cp* = --Tj nondirnens ional velocity potential • 

\|r = — V .. ^» - % nondirnens ional stream function. 

u* ~ ~- t nondirnens ional velocity in the "boundary 
0 layer. 



nondimensional velocity just outside the 
"boundary layer . 



— 2* > nondimensional "energy defect" in 
^o the "boundary layer. 

In the following computations, the asterisks are 
dropped for simplicity, and nondimensional quantities thus 
formed are always to he understood. 

Cal c ulat 1 on __o f _o_u t_e r _so 1 u t i on a - In order to calculate 
the Out er solution, it is necessary to know the variation 
of the square of the outside velocity U 2 with the veloci- 
ty potential cp along the surface. Reference 1 gives ex- 
perimentally determined values of U as a function of the 
nondimensional distance s along the surface. This rela- 
tion is plotted in figure 1. The velocity potential 

s 

. 9 = / Uds is determined by a graphical integration of 
o 

this curve. Both U and Cp are then known as functions 
of the distance along the surface. The U , <P curve, 
shown in figure 2, is det erminod by choosing corresponding 
values of U and Cp at several points along the surface. 

The next step in the calculation is the expression of 
the U ,cp curve in terms of power series. The calcula- 
tions are greatly simplified if the U 2 , Cp curve is arhi- 
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trarily divided into two regions and a separate power se- 
ries is used to approximate the curve in each region; that 

n 

is , for Cp % <p U 2 = 2 "biCp 1 

i=0 

i=0 1 

whefe Cpx is the arbitrarily chosen value of cp separating 
the two regions. It was found possible in this investiga- 
tion to obtain a satisfactory approximation to the experi- 
mental U 2 ,<p curve by using terms involving only up to 
the second power of Cp . It' appears that the U 2 , Cp curve 
should be fitted with greater care in the neighborhood of 
the separation point than elsewhere. The value of <Pi 
chosen was 0.4. The following approximate expressions for 
the U 2 ,cp curve were then obtained: 

for Cp < 0.4 U 2 = 7.380C? - 3.555C? 2 

for Cp | 0.4 U 2 = 1.393 -f 0. 414CP - 0.148Cp 2 

The approximate U , Cp curve, given by the foregoing ex- 
pressions, is also shown in figure 2. 

It is now possible to compute the energy defect zq), 
at any point in the boundary layer as a function of cp and 
if. (The subscript a) denotes the outer solution.) 

z go (»!*.) = {Po h o* * 7i<PiSi* - Y 2 cpi 2 S 2 *} 

+ i Inf** - hi*) + fcafetf * 2 72 cp 1 ( ei * + g2 *)j 

+ cp 2 |b 2 (h a - h 2 *) + p 2 h 2 * } (1) 

for cp I cp! 

y± = \ - Pi y 2 - *> 2 - P 2 

For put (3 i? (3 2 , h 0 *, h x *, h- 2 *, g x *, g 2 * 

equal to zero. The factors h„, h n , ■ h^, g, , and 
gg are junctions of y^r only, and thus do not depend on 
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the particular case. The same functions as the foregoing 

y * it ^ * 

"but formed for the variable zzz are h 0 , li 1 , h 2 

v ; ? ~ 

g x * , and g 3 * . Equation (l) corresponds to equation (19) 

in reference 2. The terras arising from the use of Dowers 
of cp greater than 2 have "been deleted. 

Substituting the values of the known coefficients in 
the expression for Zqj (equation (l)), 

z m = il,3930 h 0 * + 2.7864 g x * -I- 1.3931 g 3 *\ 

+ 9 if .3800 h x - 6.9650 h x * - 6.9556 s^..- 5.9655 g 2 *| 

+ cp 2 -1-8.8550 h 2 + 8.7070 h 2 * } (2) 

The functions g 0 a £ x , g 2 , h 0 , h x , h 2 , are computed 
from the following relations: 

So ( t) = i 0 - .(*>> i-fc 

where F(£) is the probability Integral, 

i p 2 

/ e K dp, a tabulated function* 
V TT o g 

tl fi) ^-i=u J « e (i - F (t>) 



V TT 

! 3 



*<2 



«.U) - -5 - -- 5 3 + - I 4 (i - »<{)) 

2 3*/t7 3^ 3 V 7 

h 0 U) = g 0 U) 

3a x (l) = e 0 (i) + g x (13 

h s (!) = g 0 (.;) + 2 gl (|) + g s (|) 

These fii.net ions are plotted throughout their useful rang 
in figures 3 and 4. 
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In the subsequent calculations, the power-series ex- 
pansions of these functions are required. 
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r (3) 



This computation is performed by choosing the value of Cp 
corresponding to the desired distance along the surface 
and calculating for several values of by equation 

(2). The outer limit of the boundary layer is chosen as 
the point at which relative energy deficiency — ~ C.01 

2 Z 0 

where z n = IL- . The velocity is then 0.995 the outside ve- 

Op 

locity. 

From the z m >V curves and the relations 



u = U 



the velocity profiles in the boundary layer are obtained. 
Because the integrand in the expression for y J E becomes 
infinite as the lower limit of integration is approached, 
it is necessary to evaluate the integral by a combination 
of graphic and analytic methods (reference 2, p. 10), For 



values of 



near 



that 



for small values of V 



z is expanded in a power series in \Jf , using equations 
(2) and (Z), and integrated analytically from \|/ = 0 to 
some small value of n|f , say - 0.05. From that point 
the integration may be performed graphically by :neans of 
planirneter or Coradi integraph. 

In cases whore the u,y (or z,f) curves have no 
inflection point, this procedure completes the solution 
for the particular value of cp chosen. When the u,y 
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curves have an inflection point, it is necessary to deter- 
mine certain quantities from the outer solution at the in- 
flection point to aid in fixing the boundary conditions 
for the inner solution. 



The power series in y for 2^ is again used to de- 
termine the position of the inflection point. Since the 
inflection point, in general, occurs at small values of 
\|/, powers of >Jf higher than \j/ 3 may he neglected. The 
position of the inflection point is found by equating 

This value of is denoted 



-> 2 

o z „ n 

by ijr . where the subscript j designates a joining-point 
characteristic. The values of ( 9— 



and z 1 are then 



found by substituting 'if 



GO J 



00 J 



in the power-series expansions 



for these quantities. This procedure completes the outer 
solution when the boundary- lay er velocity profiles have an 
inflection point. 

Inner so l ut ion .- The inner solution is given (refer- 
ence 2, equation (29)) in the form 




- Zr 



( 



Li 



3/4 



r 



1 - 



+ sin 



- 1 



Si 



- S 111 



+ B* 



J 



1 + B 



(4) 



where 



' 6z j. 



z = 



u 

2 



The subscript i denotes the inner solution; 
stant of integration proportional to [ B— ) 



.7=0 



I is a con- 
and z n , 



z 0 ! , and cp are derived directly from the outside pressure 
distribution, which has been experimentally determined. 
The value of £ A is taken from the outer solution and is 

equal to 



J 

- z 



CO J 
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It is to "be noted that 



S* = 0 when ^ = |j 

-at 



otner boundary condition is that g y 33 

This condition gives the equation 



9 It 
d ^> 



V 1 + B* 



An- 



(5) 



for determining B, the left-hand side of which is known 
from the outer solution. The value of B thus found is 
substituted in (4) and |j is then computed for several 

values of the ratio between 0 and 1. 

d 

It will be found, in general, that ^>±j * s n0 "^ B^qual 
to | A • In order to make the final solution continuous, 
the values of i are shifted by an amount e * * i j - 1^ j , 
that is, for the outer solution 

With this information, the £ , ijf curve can be drawn 
continuously for both solutions. As before, the u,y 
profile is determined from this curve by the following re- 
lations 



u 



V z 0 



u 



/i 



/ 



For the inner solution, this integral may be expressed 
analytically. In this case 



w n e r e 



1/4 



^2 V -z 0 ! cp v z 



o 



1 - 



- 1 



s in 



; 

i/ 1+3^ 



l n 



Fast the Joining point, the integration is most easily 
carried out graphically. 

2 

Since 1 is proportional to ] , the condition 

w / y=0 
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for separation is that 3 shall equal zero. The position 
along the surface for which this relation holds is deter- 
mined by trial and error, 

RESULTS AND DISCUSSION 



A comparison of the calculated velocity profiles and 
experimental results at several points along the surface 
of the elliptic cylinder is given in figures 5, 6, and 7. 
The calculated Pohlhausen curves, which have been discussed 
in reference 1, are shown for convenient comparison. The 
agreement becomes better forward of the separation point. 
Close to the separation point the calculations show, in 
general, too low velocities at corresponding distances 
from the surface, but the shape of the calculated separation 
point profile is in good agreement with the observed veloc- 
ity distribution at the experimentally determined separa- 
tion point. 

The positions of the calculated and experimental sepa- 
ration points are indicated on the p r e s sur e- d i s t r ibut io n 
diagrams shown in figure 8. The calculated separation 
point was found to occur at a distance along the surface 
s s= 1.92. Actual separation was observed at s = 1.99. In 
view of the many approximations and arbitrary procedures 
involved in the analysis, this agreement is considered very 
good. 

It is to bo noted that the curves given in figures 5, 
6, and 7 are independent of the Reynolds Number, except 
insofar as Reynolds Number affects the pressure distribu- 
tion. It is, of course, assumed that the boundary layer 
does not become turbulent. An example of the effect of a 
change in the pressure distribution on the position of the 
separation point is given in figure 8, where the calculat- 
ed separation point is shown for the perfect-fluid pressure 
distribution about the cylinder. In this case the position 
of the separation point is at s = 2.38, considerably aft 
the position found with the experimental pressure distri- 
bution. Figure 8 shows that, for both pressure distribu- 
tions, the separation point occurs at approximately the 
same distance, s - 0.7, aft the point of minimum pres- 
sure. Since it is to be expected that the pressure distri- 
butions for higher Reynolds Numbers will be between the ex- 
perimental curve shown and the perfect fluid-pressure dis- 
tribution, this result indicates that for all Reynolds Num- 
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bers above 24,000 the separation point for the elliptic 
cylinder will occur at a distance s == 0.7, aft the point 
of minimum pressure provided., of course, that the "boundary 
layer does not "become turbulent. 

A general view of the flow about the elliptic cylin- 
der in the l\i.A.C.A. smoke tunnel is shown in figure 9. 
The effect of separation on the entire flow configuration 
should be noted, 

CONCLUSIONS 



The computed and experimental characteristics of the 
laminar boundary layer about the elliptic cylinder are in 
good agreement. This agreement indicates that the method 
may be generally applied to the laminar boundary layer 
about any type of body provided that an experimentally de- 
termined pressure distribution is available. It appears 
that for all Reynolds Numbers above 24,000 the separation 
point for the elliptic cylinder should occur at a constant 
distance behind the point of minimum pressure, provided 
that the boundary layer does not become turbulent. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley yield, Va. , October 10, 1935. 
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Figure 2.- Experimental and approximate U^, cp curves. °^ 
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Figure 5.- Boundary- layer velocity distribution at s= 1.097. 
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Figure 7._ Boundary- layer velocity distribution near separation 
point . 




Figure 8.- Comparison of observed and calculated separation points. 
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